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ABSTRACT 

We examine the prospects for measuring the dark energy equation of state parameter w within the 
context of the still uncertain redshift evolution of galaxy cluster structure. We show that for a particular 
X-ray survey (SZE survey) the constraints on w degrade by roughly a factor of 3 (factor of 2) when 
one accounts for the possibility of non-standard cluster evolution. With followup measurements of 
a cosmology independent, mass-like quantity it is possible to measure cluster evolution, improving 
constraints on cosmological parameters (like w & 0^). We examine scenarios where 1%, 10% and 100% 
of detected clusters are followed up, showing that even a modest followup program can enhance the final 
cosmological constraints. For the case of followup measurements on 1% of the cluster sample with an 
uncertainty of 30% on individual cluster mass-like quantities, constraints on w are improved by a factor 
of 2 to 3. For the best case scenario of a zero curvature universe, these particular X-ray and SZE surveys 
can deliver uncertainties on w of ^4% to 6%. 

Subject headings: cosmic microwave background — galaxies: clusters — cosmology: theory 



1. INTRODUCTION 

Galaxy clusters have been used extensively to determine 
the cosmological matter density parameter and the ampli- 
tude of density fluctuations. Cluster surveys in the local 
universe are particularly useful for constraining a combi- 
nation of the matter density parameter f2 m and the nor- 
malization of the power spectrum of density fluctuations 
(we describe the normalization using ag, the rms fluctua- 
tions of overdensity within spheres of 8/i _1 Mpc radius; i.e. 
Henry 1997; Viana & Liddle 1999; Reiprich & Bohringer 
2002); surveys that probe the cluster population at higher 
redshift are sensitive to the growth of density fluctuations, 
allowing one to break the Qm-cs degeneracy that arises 
from local cluster abundance constraints (Eke et al. 1996; 
Bahcall & Fan 1998). Wang & Steinhardt (1998) argued 
that a measurement of the changes of cluster abundance 
with redshift would provide constraints on the dark energy 
equation of state parameter w = p/ p. 

Describing the problem in terms of cluster abundance 
only makes sense in the local universe , because, of course, 
one cannot measure the cluster abundance without know- 
ing the survey volume; the survey volume beyond z ~ 0.1 
is sensitive to cosmological parameters that affect the ex- 
pansion history of the universe- namely, the matter den- 
sity fljtf, the dark energy density VLe and the dark energy 
equation of state w. A cluster survey of a particular piece 
of the sky with appropriate followup actually delivers a 
list of clusters with mass estimates and redshifts- that is, 
the redshift distribution of galaxy clusters above some de- 
tection limit. 

Recently, it has been recognized that with current in- 
strumentation it is possible to use such surveys of galaxy 
clusters extending to redshifts z > 1 to precisely study 
the amount and nature of the dark energy (Haiman et al. 
2001). Clusters are promising tools for precision cosmolog- 



ical measurements, because they exhibit striking regularity 
and they exist throughout the epoch of dark energy dom- 
ination. Moreover, their use is complementary to studies 
of cosmic microwave background (CMB) anisotropy and 
SNe la distance measurements (Haiman et al. 2001; Levine 
et al. 2002; Hu & Kravstov 2002). Following Haiman et al. 
(2001), a series of analyses appeared that explore the theo- 
retical and observational obstacles to precise cosmological 
measurements with cluster surveys : Holder et al. (2001) 
applied the Fisher matrix formalism to the cluster survey 
problem and showed that high yield SZE cluster surveys 
can provide precise constraints on the geometry of the uni- 
verse through simultaneous measurements of He and Om- 
Weller et al. (2001) demonstrated that future SZE surveys 
might constrain the variation of the dark energy equation 
of state w(z). Hu & Kravstov (2002) examined the effects 
of cosmic variance on cluster surveys as well as including 
the effects of imprecise knowledge of a more complete list 
of cosmological parameters. Levine et al. (2002) examined 
an X-ray cluster survey, showing that a sufficiently large 
survey allows one to measure cosmological parameters and 
constrain the all-important cluster mass-observable rela- 
tion simultaneously. 

An important caveat to these works is that the authors 
assumed that the evolution of cluster structure with red- 
shift was perfectly known. In this paper, we examine the 
effects of uncertainties about cluster structural evolution 
on cosmological constraints from cluster surveys, finding 
that current survey projections that ignore this evolution 
uncertainty overstate the cosmological sensitivity of the 
survey. Furthermore, we examine the effects of survey fol- 
lowup to measure a cluster mass-like quantity M/, demon- 
strating that an appropriately designed survey can over- 
come this evolution uncertainty. 

In addition, our calculations underscore the importance 
of incorporating information from multiple observables 
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into future cluster surveys. Clusters of galaxies are dark 
matter dominated objects with baryon reservoirs in the 
form of an intracluster medum (ICM) and a galaxy popu- 
lation. Clusters can be found through the light the galaxies 
emit, the gravitational lensing distortions the cluster mass 
introduces into the morphologies of background galaxies, 
the X-rays emitted by the energetic ICM, the distortion 
that the hot ICM introduces into the cosmic microwave 
background spectrum (SZE), and the effects that the ICM 
has on jet structures associated with active galaxies in the 
cluster. These methods are largely complementary, each 
having different strengths. It appears that X-ray and SZE 
signatures of clusters are higher contrast observables than 
are weak lensing or galaxy light. That is, massive galaxy 
clusters are more prominent relative to the far more abun- 
dant lower mass halos and the large scale filaments when 
viewed with the SZE and X-ray; projection effects are a far 
more serious concern when using galaxy light or weak lens- 
ing signatures. Studies of the highest rcdshift galaxy clus- 
ters will likely be done with the SZE, because of the rcd- 
shift independence of the spectral distortion in the CMB. 
Any effort to carry out a precise cosmological study using 
galaxy clusters will undoubtedly be most effective through 
some combination of these complementary, cluster observ- 
ables. 

The paper is arranged in the following way. In §2 we 
describe two representative surveys and survey followup. 
Section 3.1 contains a description of our estimates of the 
survey sensitivity when followup is included as well as a 
description of our fiducial model. Results are presented in 
§4 and discussed further in §5. 

2. FUTURE GALAXY CLUSTER SURVEYS 

A study of using cluster surveys to probe dark energy 
begins with the redshift distribution of detectable clusters 
within a survey solid angle Aft, 

dz dzdfl ' J Q J ' dM w 

where dV/dzdfl is the comoving volume element, 
(dn/dM)dM is the comoving density of clusters of mass 
M, and f(M) is the cluster selection function for the sur- 
vey. In this analysis we take f(M) to be a step function at 
some limiting mass Mu m , which corresponds to the mass 
of a cluster that lies at the survey detection threshold. 
We use the cluster mass function dn/dM determined from 
structure formation simulations (Jenkins et al. 2001). 

In practice surveys select clusters using observables like 
the X-ray flux, SZE flux, galaxy light or weak lensing 
shear. Thus, in addition to the ingredients above, one 
requires a virial mass-observable relation (like M-L x , M— 
L sz or M-j t ). Low redshift clusters do exhibit regularity 
(i.e. David et al. 1993; Mohr & Evrard 1997), suggesting 
that observables like the ICM X-ray luminosity and tem- 
perature are good mass estimators (Finoguenov et al. 2001; 
Reiprich & Bohringcr 2002). Cluster mass to light ratios 
have been studied for decades, and it may be that this 
body of work together with modern datasets will allow 
more conclusive statements about how well galaxy light 
traces cluster halo mass. Hydrodynamical simulations lead 
us to expect that the SZE luminosity (related to the to- 
tal thermal energy within the virial region) should be the 
best ICM observable for predicting mass, but we await 



new observations with next generation SZE instruments 
to demonstrate this. 

A central feature of these mass-observable relations is 
that they evolve with redshift due to the increasing density 
of the universe at earlier times (and the changing ratio of 
distance to lense and source in the case of weak lensing). 
Within standard structure formation models, galaxy clus- 
ters form self-similarly, and so there are standard evolu- 
tion models for each mass-observable relation (e.g. Bryan 
& Norman 1998; Mohr et al. 2000; Evrard et al. 2002). 
Results to date suggest that the degree of cluster regu- 
larity locally and at intermediate redshift is comparable 
(Mushotzky & Scharf 1997; Mohr et al. 2000; Vikhlinin 
et al. 2002). However, given the central importance of 
cluster mass estimates in using surveys to study dark en- 
ergy, we can only regard these standard structure forma- 
tion models as a guide; ultimately, one needs to determine 
the evolution of cluster structure observationally. In this 
paper we examine the effects that non-standard redshift 
evolution of cluster structure would have on our ability 
to use cluster surveys to study the dark energy. As de- 
tailed below, we explore non-standard redshift evolution 
by allowing an additional dependence of (1 + z) 1 in the 
evolution of the relevant mass-observable relations. 

2.1. An X-ray and an SZE Survey 

To examine these new effects, we adopt two representa- 
tive surveys that are being promoted as ways of measuring 
the dark energy equation of state. Namely, we examine the 
following two high yield surveys: (i) a 10 4 deg 2 flux lim- 
ited X-ray survey proposed as part of the DUET mission 
to the NASA Medium-class Explorer Program, and (ii) a 
4,000 deg 2 SZE survey to be carried out with a proposed 
8m South Pole Telescope (SPT). Figure 1 contains a plot 
of the redshift distribution and limiting mass for both sur- 
veys. 

We model the DUET X-ray survey as having a bolo- 
metric flux limit of f x > 1.25 x 10~ 13 erg/s/cm 2 (corre- 
sponding to f x > 5 x 10~ 14 erg/s/cm 2 in the 0.5:2 keV 
band). For our fiducial cosmological model (see §3.2 be- 
low) this survey yields ^21,600 detected clusters, consis- 
tent with the known X-ray log A-log 5* relation for clusters 
(e.g. Gioia et al. 2001). For our mass-observable relation, 
we adopt a bolometric X-ray luminosity-mass relation 

f x {z)^d\ = A x M^E 2 (z) (1 + zf (2) 

where f x is the observed flux in units of erg s _1 cm -2 , 
dh is in units of Mpc, M200 in units of 10 15 M Q is the 
mass enclosed within a radius r2oo having a overdensity of 
200 with respect to critical and H(z) = HqE(z). where 

E 2 (z) = M (1 + zf + fl k (l + zf + fl^ 1+w) . The E(z) 
factors follow the evolution of the critical density of the 
universe p cr i t — 3H 2 /8ttG. We convert M 2 oo to M(z), the 
halo mass appropriate for our mass function at redshift z 
using a halo model (Navarro et al. 1997, hereafter, NFW; 
see discussion below regarding the effects of uncertainties 
in this conversion). Our standard evolution model ignores 
the T 1 / 2 dependence of the bolometric bremsstrahlung ra- 
diation, because X-ray surveys detect clusters using de- 
tected photons rather than detected energy. We intro- 
duce the possibility of non-standard evolution of the mass- 
observable relation with the parameter -f x . We take 7^ = 
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to be consistent with the observed weak evolution in the 
luminosity-temperature relation (Vikhlinin et al. 2002), 
and we choose (3 X — 1.807 and logA^ = —3.926, consis- 
tent with observations (Reiprich & Bohringer 2002). The 
overall h scaling of the limiting mass is ft- -1,11 . 

We model the SPT SZE survey as a flux limited survey 
with fsz > 5 mJy at 150 GHz. Within our fiducial cos- 
mological model this survey would yield ~ 13, 500 clusters 
with measured fluxes. The mass-observable relation is 



f sz (z, v)d\ = 3.781f(is)f ICM TM 200 (1 + z) 7s 



M 200 = A s 



{k B Tf 
E{z) 



(3) 



where f(v) is the frequency dependence of the SZE distor- 
tion, f sz is the observed flux in mJy, T is in Kelvin, M20Q 
is in units of 10 15 M Q , f ICM = 0.12 (e.g Mohr et al. 1999) 
and d,A is in units of Mpc (see Diego et al. 2002). We use 
log A sz = 13.466, (3 SZ = 1.48 (Finoguenov et al. 2001) and 
7sz = to model standard structure evolution. In this 
form, the overall h scaling of the limiting mass is h . 
Note that in determining the estimated uncertainties on 
cosmological parameters, we allow the normalization of 
these mass-observable relations to be free to vary. The 
survey contains enough information to solve for the best 
normalization and the cosmological parameters simultane- 
ously; therefore, shifts in model inputs like ficM within 
the observational uncertainties have minimal effect on our 
conclusions. 
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Fig. 1. — The cluster redshift distribution (heavy line) and mass 
limit (light line) of the 10 4 deg 2 DUET X-ray survey (dashed) and 
the 4,000 deg 2 South Pole Telescope (SPT) SZE survey (solid). 
The surveys are flux limited (f x > 1.25 X 10 — 13 erg/s/cm 2 and 
fsz > 5mjy at 150 GHz), and we impose a minimum cluster mass 
of 1.54 X 10 14 M Q . 

A generic problem with flux limited surveys is that at 
low redshift the implied mass limit drops well below those 
masses corresponding to galaxy clusters. The flux from 
a nearby object is spread over a much larger portion of 
the sky, and surface brightness selection effects become 
important. We model these complications by imposing a 
minimum cluster mass of 10 1 h~ l Mq. This lower limit on 
the survey mass limit is readily apparent below z ~ 0.25 
in Fig 1. 

2.2. Followup of Large Solid Angle Surveys 



The redshift distribution of clusters contains far more 
cosmological information than the surface density of clus- 
ters (Haiman et al. 2001) or the angular correlation func- 
tion (e.g. Komatsu & Seljak 2002). Thus, in both these 
surveys each detected cluster will be followed up with 
multi-band optical and near-IR photometry to provide 
photometric redshift estimates. These same data can be 
used to estimate cluster masses through their weak lensing 
effects on background galaxies (e.g. Bartelmann & Schnei- 
der 2001) and the total detected light from cluster galaxies. 

In addition, some of these clusters can be followed up 
with detailed X-ray, SZE or galaxy spectroscopic obser- 
vations that allow one to measure the mass-like quantity 
Mf(8) = M(6)/dA, which we will refer to as the followup 
mass. As an example, in the case of followup X-ray obser- 
vations that deliver the projected ICM temperature profile 
and surface brightness profile, it is straightforward to ex- 
tract the underlying ICM density p(dAd) and temperature 
profile T(dA&) to then estimate the followup mass Mf as 

k B T{6) /d\np dlnT\ 
M f{9) = ~°~q~~ (rflnf + "din?) < 4 > 

where m p is the proton mass, k B is Boltzmann con- 
stant, G is Newton constant, and the ICM number den- 
sity n = p/firrip. Note that only the shape of the ICM 
density and temperature profiles is required (i.e. knowl- 
edge of the actual distance to the system is not required). 
The followup-mass Mf can then be examined at within 
some angle 9 along with the cluster X-ray or SZE flux. 
At fixed redshift, this followup would produce an Mf-f x 
of Mf—f sz relation which would provide direct constraints 
on the structural evolution of the clusters. As we will see, 
the parameter sensitivity of these scaling relation observa- 
tions can exhibit quite different degeneracies than for the 
cluster redshift distribution, making the scaling relations 
and the cluster redshift distribution complementary. In 
§3.1 below, we describe how these survey followup obser- 
vations are included in our estimates of the cosmological 
sensitivity of the survey. 

3. COSMOLOGICAL SENSITIVITY OF A SURVEY 

3.1. Fisher Matrix Technique 

We employ the Fisher matrix technique to probe the 
relative sensitivities of two cluster surveys to different cos- 
mological and cluster structural parameters. The Fisher 
matrix information for a data set (see Tegmark et al. 1997; 

Eisenstein et al. 1998) is defined as Fij =< 



■ InC 



l 3 " dp , dp i 
dN 
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where £ is the likelihood for an observable (^p for the 
survey and Mf for the followup) and pi describes our pa- 
rameter set. The inverse F,-" 1 describes the best attainable 
covariance matrix [Cij] for measurement of the parameters 
considered. The diagonal terms in [Cij] then give the un- 
certainties on each of our parameters. In calculating these 
uncertainties, we have added the Fisher matrix for the fol- 
lowup observations (f/j ) , the Fisher matrix for the cluster 
redshift distribution (F^ ) and several external priors that 
will be discussed below. 

We construct the survey Fisher matrix Ffj following 
Holder et al. (2001) as 

ddN/dzddN/dz 1 

(5) 
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Estimated Parameter Constraints 
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where we sum over n rcdshift bins of size Az = 0.01 to 
Zmax = 3.0 and dN/dz represents the number of surveyed 
clusters in each redshift bin. The Fisher matrix for the 
followup is constructed as 



dV 

dz 
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dMf 



dn ( dM f dM f 1 



dM \ dpi dpj a 



(6) 



where Mf is a function of halo mass M and angular ra- 
dius 0, and f(dn/dM) represents the number of clusters of 
mass M for which followup mass measurements are avail- 
able in a particular redshift bin. We examine cases where 
the followup fraction / is 1%, 10% and 100%. 

To generate the followup Fisher matrix, we calculate 
the cluster binding mass within radius r — dA@ for a halo 
with virial mass M. To do this calculation we assume 
cluster mass profiles are well represented, on average, by 
NFW models with concentration index c = 5. In practice 
clusters undergo merging quite frequently and there is a 
range of halo shapes. This introduces a "theoretica" 1 un- 
certainty to the followup mass. In this analysis we take 
<Jm s — 0.3Mf to be the characteristic uncertainty in the 
followup mass measurments. This uncertainty reflects the 
observational uncertainty on individual cluster followup 
mass measurements as well as the uncertainties inherent 
in predicting the followup mass from the halo virial mass. 
As is clear from Eqn 6, the redshift and mass distribu- 
tions of the followup clusters match those of the full cluster 
survey sample; that is, we don't choose to followup only 
high redshift clusters, which would presumably provide the 
tightest constraints on our evolution parameter. We also 
choose 9 to be a dynamically varying quantity fixed to be 
95% of the virial radius corresponding to the cluster lim- 
iting mass at each redshift. Thus, followup at all redshifts 
corresponds to mass-like measurements at radii within the 
virial radius (rg(M) < r2oo(-^0)- 

3.2. Fiducial Cosmology and External Constraints 

The fiducial cosmological parameters of our model are 
h = 0.65 (i.e. Hendry et al. 2001; Ajhar et al. 2001; Reese 
et al. 2002), Vl M = 0.3 (i.e. Mohr et al. 1999; Grego et al. 
2001), Vltot = n M + n E = 1 (i.e. Netterfield et al. 2001; 
Pryke et al. 2001), w = -1, n = 0.96 (Netterfield et al. 



2001), fl B = 0.047 (Buries & Tytler 1998), and a COBE 
normalized cr 8 = 0.72 (Bunn & White 1997). Note that 
we use a rather low value of er 8 , which is consistent with 
the recent 2dF analysis (Lahav et al. 2002). Because the 
expected number density of clusters is very sensitive to 
the value of as, our fiducial SZE survey has fewer clusters 
when compared to some previous studies (i.e. Holder et al. 
2001). 

Cosmological constraints from cluster surveys are com- 
plementary to constraints from SNe la distance measure- 
ments and observations of the anisotropy of cosmic mi- 
crowave background. This is particularly true when it 
comes to using cluster surveys to measure the dark en- 
ergy equation of state parameter w (Haiman et al. 2001). 
In combination with precise CMB constraints on the cur- 
vature (fife = 0), cluster surveys enable precise measure- 
ments of the dark energy equation of state; however, when 
curvature is allowed to depart from zero- even slightly- 
the cluster constraints on w weaken considerably. For a 
flat universe , the constraints from our SZE survey assum- 
ing standard evolution and a 100% followup on w (Am) 
are 0.0406 (0.0139), whereas for = 0.01 the constraints 
are 0.1207 (0.0139). 

For the analysis presented here, we adopt relatively con- 
servative priors from future CMB anisotropy studies and 
distance measurements. We assume the power spectrum 
index n will be known to 5% , i.e. <r„ = 0.05, the Hubble 
parameter will also be known to 5%, i.e o~h — 0.0325 and 
the total density parameter fl tot will be known to 1%, i.e 
Ufc = 0.01. In addition, we take the prior on the baryon 
density parameter (Q B = 0.047) to be <tq b = 0.004. For 
reasonable values of Qb, surveys are affected by £Ib vari- 
ations only through minor effects on the transfer function 
for density perturbations (see also Levine et al. 2002). Fi- 
nally, we neglect the possibility of a variation with redshift 
in the equation of state parameter w (Weller et al. 2001). 

4. COSMOLOGICAL PARAMETER CONSTRAINTS 

Our results are listed in Table 1 and highlighted in 
the following figures. Table 1 contains a listing of la 
uncertainties on all seven cosmological and three mass- 
observable relation parameters (see Eqns 2 & 3 for def- 
initions). The first line contains a listing of the priors 
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adopted for the runs. Following that are the results for 
the SZE survey and then the X-ray survey. For each sur- 
vey we show the constraints in the case of standard evo- 
lution (i.e. 7 = 0) followed by non-standard evolution 
(7 is free parameter). The lines that follow highlight the 
impact of 1%, 10% and 100% followup. Following those 
scenarios is what we consider to be the ideal case of a flat 
universe with 100% survey followup. Finally, we show the 
case for followup only (i.e. the Fisher matrix derived from 
the redshift distribution isn't used). The constraints in 
this line provide some insights into the parameter lever- 
age that is afforded by cluster followup observations. In 
all cases the uncertainties are absolute (i.e. <tm = 0.0177 
means ft M = 0.3 ±0.0177). 
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Fig. 2. — Constraints on w and Qm for an SZE (above) and an 
X-ray survey (below). The star marks the fiducial model. Con- 
tours denote joint la constraints in five scenarios: constraints from 
dN/dz where (i— long dashed) the cluster evolution is unknown; con- 
straints from dN/dz and followup mass measurements for (ii-short 
dashed) 1% of sample, (iii— dotted) 10% of sample, and (iv-solid) 
100% of sample. The case for 100% followup plus added prior of a 
flat universe is also shown ( v-dot-dashed) . The followup mass mea- 
surements are estimated to have fractional uncertainties of 30%. 



4.1. Importance of N on- Standard Evolution 

In the standard evolution case, the SZE and X-ray sur- 
veys compare favorably, yielding la absolute uncertain- 
ties on w (n M ) of 0.1629 and 0.1659 (0.0177 and 0.0147), 
respectively. However, when one takes into account the 
possibility of non-standard evolution, the constraints on 
w weaken by almost a factor of 2 to ~ 0.25 for SZE and 
a factor 3 to ~ 0.45 for X-ray; SIm constraints weaken 
by close to a factor of 2 to ~ 0.026 for the X-ray but are 



only slightly affected in the SZE survey. The constraints 
from the cluster redshift distribution dN/dz on j sz / x are 
very weak at 0.46 and 1.29, respectively; this large un- 
certainty in the evolution of the mass-observable relation 
contributes to the weakened sensitivity to other cosmolog- 
ical parameters. 

The importance of evolution in interpreting the cluster 
redshift distribution contrasts somewhat with the results 
of the Levine et al. (2002) study, which showed that prior 
knowledge of the normalization of the mass-observable re- 
lation has only a weak effect on the cosmological sensi- 
tivity of cluster surveys (see also Diego et al. 2001). In 
their study, they only considered the standard evolution 
model. Within the context of uncertain evolution of the 
mass-observable relation one needs observations in addi- 
tion to dN / dz to determine the evolution parameter 7 and 
regain high sensitivity to the equation of state parameter 
w. Next we examine the effects of including followup mass 
measurements . 

4.2. Effects of Followup Mass Measurements 

Figure 2 contains joint constraints on Q m and w for the 
two surveys that highlight the effect of survey followup. 
For each survey we show constraints with non-standard 
evolution and no followup (long dashed), 1% (dashed), 
10% (dotted), and 100% (solid) followup along with 100% 
followup in a flat universe (dot-dashed) . The figure makes 
clear that even a limited followup program can greatly im- 
prove cosmological constraints. 

Table 1 shows that as the followup fraction increases, the 
constraints on w in the SZE survey tighten from 0.25 (no 
followup) to 0.14 (1% followup) to 0.12 (10% and 100% fol- 
lowup). The difference between 10% and 100% followup 
is minimal, suggesting that the ten times higher cost of 
full survey followup is not a worthwhile investment when 
viewed solely from the perspective of obtaining constraints 
on the equation of state of the dark energy. 

Notice that the constraints on w in the cases of even 
limited followup are somewhat better than the constraints 
in the cases where we assume complete knowledge of the 
evolution of the mass-observable relation. In the X-ray 
survey, a program to followup as few as 1% of the clus- 
ters can offset the increase in uncertainties that we see in 
going from standard evolution to non-standard evolution. 
In the SZE survey, 1% followup produces constraints that 
are somewhat better, reducing the uncertainty on w from 
0.1629 in the standard evolution scenario to 0.1404 in the 
non-standard evolution + l%followup. This can be traced 
to our assumption in these calculations that the redshift 
distribution of followup mass measurements matches the 
redshift distribution of the full survey. The higher redshift 
followup measurements contain more information about 
evolution, and the SZE survey probes to higher redshift 
than does the X-ray survey (see Fig. 1) 

Table 1 contains a listing of the effects of followup on 
all parameters. It is clear that followup mass measure- 
ments dramatically reduce the projected uncertainties on 
cosmological and scaling relation parameters. As is ev- 
ident from the last column in the table, even a modest 
followup of 10% of the clusters reduces the uncertainty on 
7 from 0.46 to 0.10 for the X-Ray case and 1.29 to 0.12 
for the SZE survey. With full followup, the constraint on 
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w shrinks from 0.25 to 0.12 in the SZE and 0.45 to 0.14 in 
the X-ray survey. Even with followup of only 1% of the 
clusters in the SZE survey, one reduces the uncertainty on 
w by half. Comparison of the uncertainties on 7 for non- 
standard evolution with no followup to the followup only 
case shows clearly that followup is very effective in con- 
straining 7. For example, for the X-Ray case the followup 
itself can constrain 7 to ~ 0.19 compared to 1.29 that one 
can achieve from the survey only. This underscores the ad- 
vantage of having a followup program for cluster surveys. 
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Constraints on the mass-obscrvable relation normal- 
ization A and redshift evolution (1 + z) 7 for an SZE (above) and 
an X-ray survey (below). The star marks the fiducial model. Con- 
tours denote joint la constraints in four scenarios: constraints from 
dN/dz where (i— long dashed) the non standard cluster evolution is 
unknown; constraints from dN/dz and followup mass measurements 
for (ii— short dashed) 1% of sample, (iii-dotted) 10% of sample, and 
(iv— solid) 100% of sample. The followup mass measurements are 
estimated to have fractional uncertainties of 30%. 

In Fig. 3 we show the constraints on the mass-observable 
relation normalization A and the evolution parameter 7 for 
the four cases: no followup (long-dashed) and a followup 
of 1% (short-dashed), 10% (dotted) and 100% (solid) of 
the clusters. Followup has strikingly different effects in 
the SZE and X-ray surveys. Followup in the SZE survey 
is more effective at constraining the evolution parameter 
7 SZ , due to the greater redshift depth of this survey. The 
differences in the constraints on log 10 A generally reflect 
the different definitions of the normalization and its rela- 
tionship to halo mass (see Eqns 2 & 3). In contrast to 
the previous figure that shows the effects of followup on 
w constraints, it is clear from this figure that if one really 
wants to understand the normalization and evolution of 
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Fig. 4. — The sensitivity of the Followup (above) and Survey (be- 
low) to change in ui is shown for both X-Ray and SZE surveys. The 
four cases are: (i-solid) SZE survey weighted sensitivity of Fij to 
change in w and (ii-dashed) X-Ray survey weighted sensitivity Fij 
to change in w (see Eqn 5 & Eqn 6 for an explanation of Fij); (iii- 
dot-dashed) same as (i) but per unit cluster and (iv-dotted) same 
as (ii) but per unit cluster. The redshift distribution of the clusters 
are shown in. 

4.3. Redshift Variation of Parameter Sensitivity 

In Fig 4 we display an estimate of the redshift variation 
of the survey and followup sensitivity to w. We do this 
by examining the derivative with respect to redshift of the 
w-w components of the survey and followup Fisher matri- 
ces. These are shown using heavy lines for both the X-ray 
(dashed) and SZE (solid) surveys. We also show an esti- 
mate of the per unit cluster sensitivity using lighter lines 
for the X-ray (dot-dashed) and SZE (dotted) surveys. The 
axes for the heavy curves are included on the left of the 
figure, and the axes corresponding to the per unit clus- 
ter values are located on the right of the figure. In both 
the X-ray and SZE cases and for both followup and the 
survey itself, the heavy curves show sensitivity that peaks 
at lower redshift than in the lighter curves that show the 
per unit cluster sensitivity. This is simply a reflection of 
the redshift distributions of the clusters detected in the 
surveys, which peak at z < 0.5 in both surveys. 

Consider now the lower panel. The higher redshift na- 
ture of the SPT SZE survey relative to the DUET X-ray 
survey is clear in this panel, where we see that the w sen- 
sitivity of the SZE survey peaks near z ~ 1 (solid line). 
For the X-Ray survey, the w sensitivity peaks at z ~ 0.7 
(dashed line). However, it is clear that the sensitivity per 
unit cluster ((dF^ w j ' dz) j (dN j dz) , denoted by dotted line 
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for SZE and dot-dashed line for X-Ray) increases as we 
go to higher redshift. This emphasizes the importance of 
high redshift cluster surveys for probing w. High redshift 
clusters in the DUET survey are even more sensitive to 
w, because they are more massive and lie well beyond the 
exponential cutoff in the mass function; however, these 
clusters are so rare that none are detected in the DUET 
survey. 

The sharp cutoff and first mini-peak at redshift z ~ 0.3 
of dF^ w /dz for the survey is a direct result of our re- 
quirement that clusters have masses above 10 14 h~ 1 M Q . 
In Fig 1 the plot of the limiting mass becomes flat at red- 
shifts below z ~ 0.3; we include this mass cutoff for several 
reasons: (1) the mass-observable scaling relations we've 
adopted are for clusters, and they are inappropriate for 
group scale systems and (2) a flux limited survey becomes 
sensitive to surface brightness limitations at low redshift, 
where the flux from these nearby objects is spread over a 
larger and larger solid angle. Introducing a minimum mass 
in our survey causes this interesting artifact in the to sen- 
sitivity at low redshift, which can be understood by con- 
sidering the competitive behaviour of the w-dependences 
of the limiting mass, the survey volume clement, and the 
growth factor of density perturbations. As long as the 
limiting mass is constant, the opposite sensitivities of the 
volume element and the growth factor to w determine the 
net to-sensitivity of the survey. At low redshift, the w- 
sensitivity of the volume element dominates over that of 
the growth factor. Beyond redshift ~ 0.3, the limiting 
mass suddenly rises above the minimum, allowing the net 
w sensitivity to include that of the limiting mass. The lim- 
iting mass is sensitive to w primarily through the angular 
diameter or luminosity distance to the redshift of inter- 
est, and this sensitivity combines with that of the growth 
factor to offset to a larger degree the w sensitivity of the 
volume element. Note that dF^ w /dz is positive definite 
which is the reason for the visual appearance of a break in 
the sensitivity, which is actually reflecting an underlying 
change in sign of the to sensitivity. 

For the followup (upper panel), the sensitivity per unit 
cluster peaks at z ~ 0.8 for both the surveys. However, 
survey weighted sensitivities of the followup to variation 
in w are peaked at much lower redshifts, where the sur- 
vey yields are much higher. This is due to the fact that 
the to-sensitivity of the followup is a balance between the 
to-sensitivities of the mass observable and the number of 
clusters one has at that redshift to make the measure- 
ment. The fact that the redshift dependence of the cluster 
mass observable is similar for both X-ray and SZE just 
reflects the weak mass dependence of this sensitivity. Sur- 
vey strategists should consider a followup program that 
targets predominantly redshift z ~ 1 clusters if their goal 
is to constrain the evolution parameter and improve con- 
straints on to; however, having evolution information over 
the entire redshift range of the survey is critical to testing 
the form of the non-standard evolution model, which we 
parametrize here simply as (1 + z) 1 . 

5. DISCUSSION AND CONCLUSIONS 

Any attempt to precisely measure the dark energy equa- 
tion of state w with cluster surveys will require (i) a strong 
external prior on the curvature (presumably from CMB 
anisotropy studies) and (ii) an understanding of the evolu- 



tion of the relation between cluster halo mass and observ- 
able properties like the X-ray luminosity, SZE luminos- 
ity, galaxy light or weak lensing shear. We have examined 
the effects of current uncertainties about cluster structural 
evolution; for two recently proposed cluster surveys the 
estimated constraints on to are ~2-3 times weaker than 
if one assumes full knowledge of cluster evolution. Con- 
straints on other interesting cosmological parameters are 
also weakened (see Table I). 

Followup observations to measure cluster masses di- 
rectly will enable one to solve for cluster structure evo- 
lution and to enhance cosmological constraints. We have 
examined the effects of followup mass observations from 
hydrostatic or dynamical methods, and we find that even 
modest followup of 1% of the cluster sample can improve 
survey constraints. Full followup with mass measurements 
that are 30% uncertain, on average, provide cosmological 
constraints that match or surpass those possible through 
dN/dz alone with full knowledge of cluster evolution. Full 
followup with weak lensing mass measurements is cur- 
rently being planned for the SPT SZE survey. 

The implications are quite interesting. Essentially, to do 
precision cosmology with cluster surveys and followup, we 
need only know that clusters conform to mass-observable 
scaling relations, and that these relations evolve in some 
well behaved manner. Then, together with our well estab- 
lished theoretical framework for structure formation, the 
observed cluster redshift distribution and followup masses 
of as few as 1% of the sample then provide enough in- 
formation to deliver precise constraints on cosmological 
parameters and the character and evolution of the mass- 
observable relation simultaneously. In a sense cluster sur- 
veys with limited followup arc self calibrating: one gains 
detailed knowledge of the structure of the tracers (i.e. 
galaxy clusters) and detailed knowledge of the evolution 
of the universe from the same dataset. 

We have focused here on the mean equation of state pa- 
rameter to, and for the two surveys considered, we examine 
the redshift variation of the sensitivity to to. In the SPT 
SZE survey, the sensitivity to the dark energy equation of 
state peaks at z ~ 1 with sensitivity at or above half the 
peak for 0.65 < z < 1.5. In the DUET X-ray survey, the 
sensitivity peaks at z ~ 0.7 with sensitivity at or above 
half the peak for 0.45 < z < 1.0. In the case of both 
surveys the to sensitivity of followup mass measurements 
peaks near redshift z ~ 0.35. In the case of the clus- 
ter redshift distribution, the most information about to is 
provided by the highest redshift clusters, and so deeper, 
more sensitive surveys will in general be better for studies 
of the dark energy equation of state. 

One interesting feature of our analysis is the orientation 
of the elliptical constraints on to and Qm (see Fig. 2). In 
general, the rotation of the parameter degeneracy can be 
understood as the result of competing effects of changes in 
the volume element and the growth factor as parameters 
vary. Variations in to (and SIm) affect the survey yield 
in different ways at different redshifts, and so the w-CIm 
degeneracy depends on the redshift distribution of a par- 
ticular survey. Rotations of parameter degeneracies occur 
as the maximum redshift of the survey is varied (Levine 
et al. 2002; Hu & Kravstov 2002). We have further found 
that changing the prior on Q tot and changing the degree 
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of mass followup on a survey also result in rotations of 
the parameter degenaracy. This behavior has interesting 
implications for the design of cluster surveys that are op- 
timally complementary to CMB anisotropy and SNe la 
distance measurements, and it deserves further study. 

In addition, we emphasize that the final constraint 
on the determination of cosmological parameters depends 
sensitively on the survey strategy and also the details of 
the followup. For example, a different definition of Mf(6) 
would lead to slightly different uncertainties. Changing 
9 from a quantity that varies with redshift to some fixed 



value leads to modest variations of the constraints. In gen- 
eral, best results can be obtained by optimizing Mf(9) so 
that one probes as mush of the virial mass as possible. 
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